Self-Energy Correction to the Hyperfine Splitting for Excited States 
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The self-energy corrections to the hyperfine sphtting is evaluated for higher excited states in 
hydrogenlike ions, using an expansion in the binding parameter Za, where Z is the nuclear charge 
number, and a is the fine-structure constant. We present analytic results for D, F and G states, 
and for a number of highly excited Rydberg states with principal quantum numbers in the range 
13 < n < 16, and orbital angular momenta £ = n — 2 and I — n — 1. A closed-form, analytic 
expression is derived for the contribution of high-energy photons, valid for any state with £ > 2 and 
arbitrary n, £ and total angular momentum j. The low-energy contributions are written in the form 
of generalized Bethe logarithms and evaluated for selected states. 
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I. INTRODUCTION 



The self-energy correction to the hyperfine splitting 
is the dominant quantum electrodynamic (QED) correc- 
tion to the magnetic interaction of the bound electron 
with the field of the nucleus. The hyperfine interac- 
tion energy of electron and nucleus is proportional to 
gN a{Za)'^m'^/mi^ , where gN is the nuclear g factor, and 
TOe and totv arc the electron and nuclear masses, re- 
spectively. Rclativistic corrections enter at relative or- 
der (Za)^. The dominant QED correction is due to the 
anomalous magnetic moment of the electron and enters 
at relative order a. Here, we consider the QED correction 
of order a {Za)^ , which is the sum of a high- and a low- 
energy part. Relativistic corrections to the anomalous 
magnetic interaction give one of the dominant contribu- 
tions to the high-energy part, which can otherwise be 
calculated on the basis of a form-factor approach, using 
a generahzed Dirac equation in which the radiative effects 
and the hyperfine interaction are inserted "by hand." 
The low-energy part constitutes a correction to the Bethe 
logarithm due to the hyperfine interaction. It can be for- 
mulated as a hyperfine correction to the self-energy, the 
effect being equivalent to the self-energy correction to the 
hyperfine splitting mediated by low-energy virtual pho- 
tons [up to order a [ZaY]. 

In our treatment, we follow the formalism of nonrel- 
ativistic QED (NRQED) detailed in Ref. Q, and refer 
to Refs. 0-01 for a number of previous investigations re- 
garding the treatment of the self-energy correction to the 
hyperfine splitting in systems with low nuclear charge 
number. 

Our paper is organized as follows. The general formal- 
ism of the hyperfine interaction is described in Sec. |lll 
For the self-energy correction, the low-energy part is 
treated in Sec. IIIIl and the high-energy part is calcu- 
lated in Sec. IIVI Results and theoretical predictions are 
discussed in Sec.|Vl Conclusions are reserved for Sec. lVIl 
Natural units with ft = c = eg = 1 arc used throughout 
the paper. 



II. FORMALISM 

Following the derivation in Ref. Q, the magnetic 
dipole field of the nucleus is described by the vector po- 
tential 
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where x is the coordinate vector and r = \x\. The curl of 
this vector potential yields the magnetic field 



Ms = V X Ams{x) - -- pdix) - ^^^3 ^ , (2) 



and the fully relativistic hyperfine interaction Hamilto- 
nian thus reads 

^ ? e ^ fix X e ^ X X d 

Hms = ~ea • AMs[xj = —a — = — /i — (3) 

47r r-^ 47r r-^ 

The hyperfine interaction couples Dirac eigenstates to 
the magnetic field of the nucleus. The electronic states 
can be written as \njjn) = \nj£m), where n is the prin- 
cipal quantum number, and the orbital and total an- 
gular momenta of the electron {£ and j, respectively) 
can be mapped to the Dirac angular quantum number 

K = {—iy^^^2 (j + i). Finally, m is the projection of the 
total electron angular momentum onto the quantization 
axis. In this article, we sometimes suppress the orbital 
angular momentum £ in the notation because we consider 
the coupling of the total electron angular momentum j 
to the nuclear spin. Nuclear states are denoted as \IM), 
where / is the nuclear spin and M its projection onto 
the quantization axis. They are coupled to the electron 
eigenstates Injm) by the hyperfine interaction, to form 
states with quantum number \nfmflj) which are eigen- 
states of the total Dirac -l-hyperfine Hamiltonian (/ is the 
total electron-l- nuclear angular momentum, and m/ is its 

projection). Using Clebsch-Gordan coefficients Cj^/-^, 
the \nfmflj) states can be written as 



Wfruflj) 
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The hyperfine energy Ai?hfs thus reads 

£^iifs = {nfmflj \Hms\ nfmflj) 



(5) 



Using the Wigner-Eckhart theorem, the hyperfine energy 
can be rewritten as 

Ems = [/(/ + 1) - /(/ + 1) - jij + 1)] 



\x X a L 
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where is the Dirac eigenstate with a definite an- 

gular momentum projection and [x x d]g is the z 
component (zero component in the spherical basis) of the 
indicated vector product. 

We have thus separated the nuclear from the electronic 
variables. A detailed analysis of the separation of the 
nuclear variables can also be found in Ref. 0. This 
procedure allows to reduce the evaluations of the hyper- 
fine structure and corrections to it, to the evaluation of 
matrix elements of operators acting solely on electronic 
states. Specifically, we consider corrections to the state- 
dependent electronic matrix element Gg, where 



TOe?"' 



The hyperfine interaction energy thus is 

Ems = [/(/ + !) - /(/+!) - + 

I rriN 



(7) 



(8) 



Rclativistic atomic theory leads to the following result 
for Oe (see Refs. 0,11) 

n (7 ^3 n{2K{^ + n-\K\)^N) 



where 7 = -^/k^ — (Za)'^. The effective principal quan- 
tum number is iV = {n — \k\)^ + 2{n — |k|)7 + k^. 



III. LOW-ENERGY PART 

A. General Formalism 

In order to treat low-energy virtual photons, we apply 
a Foldy-Wouthuysen transformation to the total Hamil- 
tonian Ht which is the sum of the Dirac-Coulomb Hamil- 
tonian and the relativistic hyperfine interaction Hamilto- 
nian. 



Ht = Hd+Hms = d-p + prrie 



Za 



'ed-Ahfe(f). (10) 



The Foldy-Wouthuysen transformation of this Hamilto- 
nian is carried out as described in Refs. [tIIsI fl0l - [T2j . The 



only difference to the case of the ordinary Dirac Hamilto- 
nian is that the odd operator O used in the construction 
of the transformation [ll| now reads 

= d-p~ea-AMsix) (11) 
instead of d? • p. The result of the transformation. 



H' + UHtU- 



FW 



HFS : 



(12) 



is the sum of the Foldy-Wouthuysen Hamiltonian ffpw 
from Ref. [lol | , and i^HPs is the nonrelativistic hyperfine 
splitting Hamiltonian [2, 

It consists of the three parts, 
r 47r 



(7(5 (x) . 
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(14a) 
(14b) 
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whose designation is inspired by the apparent angular 
momentum association of the terms. Following the nota- 
tion in Ref. Q , lowercase letters in the subscript are used 
to label relativistic operators, whereas nonrelativistic op- 
erators are denoted by uppercase letters in the subscript. 
However, we use the lowercase notation for the scaled 
vector quantity h in order to denote the electronic oper- 
ators in the nonrelativistic hyperfine Hamiltonian. Fur- 
thermore, X = is the position unit vector. The zero 
component (z component) /iq = ^S,o + ^c.o + ^i.o of the 
Hamiltonian "vector" h therefore reads as 



u ^'^ xr^\ , 3(0" ■S:)xo- (To 
^0 = iT—^<^od{x) + 



3to2 



(15) 



With the help of Iiq, the nonrelativistic limit of Eq. 
is obtained as 



\k\ n^{2K + l){K^ ~ i) 
which we use to define the nonrelativistic quantity 

NR 9N 
Ef — ^HFS — EhFS ~ a— 



(16) 



(17) 



NR 



X [/(/ + 1)- J(/+l)-j-(j + l)]0 

which is commonly referred to as the Fermi energy. The 
relativistic and QED corrections can be expressed as mul- 
tiplicative corrections of 0e, via the replacement 



qNR ^ qNR ^ 50rcl ^ ^qQED] 



(18) 
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By expanding to second order in Za, we obtain 



6ef = {Za) 



I2n^ - 1 3 1 

2k2(2k- 1)(2k + 1) ^ 2^H 

3 -8k 



2n2(2K- 1)^ 
and the corresponding energy shift 



The QED term 
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The hyperfine correction to the interaction current is 
then given as 







n' j' I'm' 
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{n'j'£'m'\Sjl^Ps\nj£^) 



iiZaf ^ 

n' j' I' 7n' 



me(ZQ;)2 
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is the subject of this paper. For the QED corrections 
terms up to relative order a{Zaf' with respect to the 
nonrelativistic hyperfine splitting will be considered. In 
order to do the calculation, we need the three terms from 
Eq. ([T^ and a further correction to the electron's tran- 
sition current, due to the hyperfine interaction. Namely, 
in the presence of the hyperfine interaction, the kinetic 
momentum of the electron finds a modification 



P P A P 
> ^hfs = 

m„ rUp m„ nip 



The current 

has the zero component 



|e| 
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{-yex + xiy) , 
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The term containing the logarithm of e, which is a scale- 
separation parameter that cancels when high- and low- 
energy parts are added [l3j . vanishes after angular inte- 
gration in the matrix element. The structure of the loga- 
rithmic term here is very similar to the Bethe logarithm 
encountered in Ref. [lj|- Terms of this form will arise 
for the other corrections in the low-energy part as well. 
In the following, these terms are denoted as /3hfs and 
are evaluated numerically with the methods described in 
Ref. [l^. Thus, the low-energy correction due to the 
nuclear-spin dependent current is 



HFS ■ 



(27) 



Next, we treat the corrections to the Hamiltonian, to 
the energy and to the wave function. The perturba- 
tion due to the hyperfine splitting Hamiltonian yields the 
term [we define the resolvent G{uj^) = 1/ {En — Hs — i^^:)] 



which is used in the calculations below. 



B. Specific Terms 

Following Ref. Q, there are four corrections, which 
arise from the correction of the interaction current, from 
the correction of the Hamiltonian, from the correction of 
the reference state energy, and finally from the correction 
of the reference-state wave function. We first treat the 
hyperfine correction to the interaction current and to this 
end, define a useful normalization factor 



1 



1 



(25) 
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+ -{Zari3'^%, 
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The correction to the energy denominator in the 
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Schrodinger propagator can be written as 



2aN 
37r 



2 p 



nj4 



In 



X {njq \ho\ njq) + ^{Zafpfis , (29) 



where the prime indicates the reduced Green function. 
Finally, the correction to the wave function due to the 
hyperfine splitting Haniiltonian is 

4a7V 



Sir 



^k 




P r<( \ P 



1 



En — Hs 



(30) 



Using commutator relations, one can finally sum up all 
four corrections in the low-energy part to 



(59, 



a 



a 



In 



TT 



(Za)2/3] 



TOe(.^a)^ 
HFS , 



{njl^\[p\[h^,f]]\njl\) 



where /3hfs is the sum 



Phfs = Phfs + /'HFS + Phfs + Phfs ■ 



(31) 
(32) 



The double commutator 



5Qi 



njq) (33) 



vanishes for states with £ > 2 up to and including order 
(Za)^, and hence SQl takes the very simple form 



-(Zq;)^/3hfS 



(34) 



IV. HIGH-ENERGY PART 



Up to relative order a{Za)^Ep, it is sufficient 



to 



consider the problem on the level of the modified Dirac 
Hamiltonian 



H 



(m) 



a 



p-eFi{y^)A +/3me + Fi(V2)T/ 



2me V 



[i^-E-/3^-B , 



(35) 



where i^i and F2 are the one-loop Dirac and Pauli form 
factors of the electron, respectively. Their expressions 
are known (see Chapter 7 of Ref |16j). 

The Fi form factor slope gives rise to the following 
effective interaction 



-eFi{0)V^d-AMs 



3n 



In 



2e 



11 
24 



V^Fhfe. (36) 



Up to the order a{Za)^ Ep, we may write the correction 
in terms of the nonrelativistic hyperfine Hamiltonian ho, 



SQ 



H.l 



(37) 



However, as already pointed out, the matrix element of 
V^/iQ vanishes on states with £ > 2 which are relevant to 
our investigations, and so 



5Q 



HA 



0. 



(38) 



The second correction is a second-order perturbation in- 
volving the Fi correction to the Coulomb potential. 



H.2 




Again, V^V is proportional to the Dirac 5 and therefore 
vanishes for states with £ > 1. Accordingly, for states 
with ^ > 2 we have 



5Q 



H,2 



0. 



(40) 



The Pauli F2 form factor gives rise to a second-order 
perturbation involving a magnetic moment correction to 
the Coulomb potential. 



dQH,3 = 2AAF2(0) 



X ( 



(41) 



2m, 



■7- vy 



E,h — Ht 



hfs 



where F2 is the magnetic form factor. For ^2(0), the 
Schwinger value ^2(0) = ^ may be used. After a Foldy- 
Wouthuysen transformation, we can write 5Qh,3 as the 
sum of two terms. The first of these, 6QH,3n, involves no 
mixing of upper and lower components in the Dirac wave 
function and reads 



(56 



ZOL 



En — Hf. 



HFS 



njq ) . (42) 
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We find the following general result for states with ^ > 2, 



a 



( 1 60£^ + 120^3 _^ 55£2 _ 5^ _ 3 



1 



2nAt(2£ + l) 8n2K(£+|)(^-i) 



(43) 



Lower components of the Dirac wave function give rise 
to the mixing term 



aJ\f —K 



2m„ 



We find the general result 



nj4 



nji\ ) . (44) 



AO - ^2 1^1 (2«+l)(/.^-|) 

^e,,3™ - -(Za) 4^.(^. + i) (2,-l)(2^ + 3)(£+i) 



1 



n2 £(€+1) 



(45) 



The F2 correction to the magnetic photon exchange of 
electron and nucleus gives rise to the effective interaction 



-^2(V^)— /3S-Bhfs = F2(V^) —±l3fl-{hs + ha 



Here, and are the generalizations of hs and ho to 
4x4 matrices, 



47r ^ 
hs = S (5(f) , 

j> 3 (S • a;) z — S 
Z mi r'^ 



(47a) 
(47b) 



Taking F2{V'^) w ^2(0) in Eq. (gll), we obtain the cor- 
rection 

dQn.i = AAF2(0) {^P 1/3 (/i«,o + /^d,o)| ^) 

= ^ {'^\P{hs,o + hd,o)\^) . (48) 

Generalizing results from Ref. for the term of relative 
order a, we find the result 



(56 



HA 



J_ , / 24^^3 + 18^.^-^.-1 

4k ^ V8«;3(4k3 + 4k2 -K- 1) 
3 1 1 1 1 - 3k 



8n \k\k -n? 2k 2k — 1 
Taking the slope of F2 in Eq. pS)) . we obtain 

^^2(0);7^/3V2S-Sms 
2mp 



(49) 



12^ 



eTOc 
47r 



(50) 



TABLE I: Low-energy contribution /?hfs of the self-energy 
correction for the hyperfine splitting for D states = 2). 



/Shfs {nD; 



3/2j 



/3HFs(n-D5/2) 



-2.068 39(5) X 10"^ 
-1.302 99(5) X 10^2 
-1.025 92(5) X 10"^ 
-0.943 28(5) X 10"^ 



-3.455 22(5) x 10"^ 
-3.793 94(5) x 10"^ 
-3.965 95(5) x 10"^ 
-4.084 82(5) X 10"^ 



TABLE 11: Low-energy contribution /3hfs of the self-energy 
correction for the hyperfine splitting for F states = 3). 



^HFs(nf5/2) 



/3hfs [nF-j/2 



-1.021 46(5) X 10"^ 
-0.683 94(5) X 10"^ 
-0.504 64(5) X 10^2 
-0.407 32(5) X 10^2 



-0.953 04(5) X 10-2 
-1.077 62(5) X 10~2 
-1.141 28(5) X 10-2 
-1.182 43(5) X 10-2 



TABLE IIL Low-energy contribution /?hfs of the self-energy 
correction for the hyperfine splitting for G states = 4). 



^HFs(nG7/2) 



/3hfs(?iG9/2) 



-0.368 23(5) X IQ-^ 
-0.147 39(5) X 10-2 
-0.012 93(5) X 10-2 
0.072 16(5) X 10-2 



-0.341 76(5) X 10-2 
-0.388 58(5) X IO-2 
-0.414 24(5) X 10-2 
-0.430 98(5) X 10-2 



with F^(0) = Q;/127r. As F^(0)V2 already is of relative 
order a[Za)^ , this operator only has to be applied to the 
nonrelativistic wave function where it vanishes for states 
with t> 2 and thus 



sen, = ^ U 



{hs.Q + hofl) 



0. 



Finally, since 

,56^.1 = (5eH,2 = (5eH,5 = , 
we have for the high-energy part the result 

5Qh = 5QH,3n + SQH,3m + SQhA 



(51) 
(52) 

(53) 



It is quite surprising that the result obtained by adding 
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TABLE IV: Low-energy contribution /?hfs of the self-energy correction for the hyperfine splitting 
for highly excited states. The numbers in parentheses are standard uncertainties in the last figure. 
The exponent of the numerical data is chosen to be the same as in Tables HI nil and|TITl illustrating 
the monotonic decrease of the coeffcients with the angular momentum quantum numbers. 



n 


I 


2i 


K 


/3hfs 




2i 




/3hfs 




16 


15 


29 


15 


0.006 310(5) X 10" 


-2 


31 


-16 


-0.002 130(5) X 10" 


-2 


16 


14 


27 


14 


0.016 397(5) X 10" 


-2 


29 


-15 


-0.003 041(5) X 10" 


-2 


15 


14 


27 


14 


0.006 888(5) X 10" 


-2 


29 


-15 


-0.002 795(5) X 10" 


-2 


15 


13 


25 


13 


0.019 372(5) X 10" 


-2 


27 


-14 


-0.004 086(5) X 10" 


-2 


14 


13 


25 


13 


0.007 420(5) X 10" 


-2 


27 


-14 


-0.003 741(5) X 10" 


-2 


14 


12 


23 


12 


0.023 029(5) X 10" 


-2 


25 


-13 


-0.005 617(5) X 10" 


-2 


13 


12 


23 


12 


0.007 794(5) X 10" 


-2 


25 


-13 


-0.005 122(5) X 10" 


-2 



the above expressions, 
a I 1 



5Q 



H 



TT An 



+ {ZaY 



1 24k3 



18^2 



1 -f 120£3 -f- 55£2 -5^-3 

ni^ {21+ If (4£3 + 8£2 £ - 3) 



J + 



(2k-H 1)(k 



2 _ 1^ 

4^ 



'+l)4i(j + l) (2£-l)(2£ + 3)(£+i) 



1 3 4(j + l) + (2£+l) 
n 8k (2£+l)(j + i) 
(2j + l)(2K- 1)(2k-H1)^ 



1 

8^ 



1 — 2k k 
3£(£+1) 



2j{j + 1)(2£ + 3) (4£2 - 1) +!)(€- i) 



can actually be simplified quite considerably, 



a 



1 



+ 



TT [4k 

1 4k- 1 



{Za) 



13k 6k + 1 



n 8 |k| K2(2K-hl) 
(4k-K 1)(6k-H 1)(6k2 -h3K- 1) 
n2 2k(1 - 2k) ^ 8k3(2k + 1)2(2k- l)(K-hl) 



(54) 
(55) 



V. RESULTS AND PREDICTIONS 

The total self-energy correction to the hyperfine split- 
ting is obtained as the sum of the high- and low-energy 
parts given in Eqs. ((55|) and ((34)) . which reads 



5Q 



+ {Za) 



a \ I 

TT \4k 

1 (4k+1)(6k+1)(6k2 + 3k-1) 
8^ (2k + 1)2(2k - 1)(k + 1) 



(56) 



1 4k- 1 



13k 6k + 1 
'n Syi^ k2(2k-1- 1) ^ ^2k(1 - 2k) 



HPS 



The Bethe logaritm type correction /3hfs is given in 
Eq. dsn). 



Restoring the reduced-mass dependence [we define 
r{M) = TTie/mTv] and adding the relativistic correction 
of relative order (Za)^, we find that 



RqqC 



9N 



[l + r(A/-)]3|K|(2K+l)(K2-i) 

X [/(/ + 1) -/(/+!)- J (.7 + 1)] X |l + (Za)2 
12k2 - 1 3 1 3 - 8k 1 



2k2(2k- 1)(2k+1) 

a 1 a , , T 3 k 
+ -— + -{Za) 

TT 4k TT 



3 1 

2n |k| 
6k + 1 



2n2(2K- 1) 



n |k| k2(2k + 1) 
1 (4K-H)(6K-t-l)(6K^-F3K- 1) 
8k3 (2k+ 1)2(2k- 1)(k + 1) 



1 4k - 1 
^2k(1-2k) 

' /^HFS I ■ 

(57) 



Numerical data for /3hfs for D, F, and G states, and 
selected Rydberg states, can be found in Tables HI HIllIIIl 
and Table IIVI respectively. The latter are relevant 
for a possible determination of fundamental constants 
from Rydberg state spectroscopy in hydrogenlike ions of 
medium nuclear charge number (see Ref . [18| ) . 

A numerical example: For the transition |1) o |2) in 
atomic hydrogen {Z = 1) where 

|1) = |n = 15,£= 14,j = 29/2,/= 15) , (58) 
|2) = In = 16, £ = 15, J = 31/2, / = 16) , (59) 

we find the following frequency shift from the Dirac value 
indicated in Eq. (2) of Ref. 

Ai/hfs 1^2 = 96.759 8630(8) Hz + 78.764 693(13) Hz , 

(60) 

where the first term is due to the hyperfine effects cal- 
culated here, and the second term is due to relativistic 
recoil and QED effects calculated in Ref. [l^- The final 
theoretical prediction for the shift from the Dirac value 



7 



is 

Aj^hfs,i^2 = 175.524 556(13) Hz , (61) 

where the fundamental constants of CODATA 2006 [Hi 
have been used in the numerical evaluation. The next 
higher-order term neglected here is the recoil correction 
of relative order {Zol)^ ^{M), for which a general expres- 
sion has been derived in Ref. (20j [the corresponding ex- 
pression also is given in Eq. (42) of Ref. Q]. The recoil 
correction is numerically suppressed for Z = \. 

VI. CONCLUSIONS 

Rydberg states of hydrogenlike ions with medium nu- 
clear charge number have been proposed as a device for 
the determination of fundamental constants [l^. Here, 
we demonstrate that it is possible to obtain accurate 
theoretical predictions for transition frequencies even in 
cases where the nucleus carries spin. To this end, we cal- 
culate the self-energy correction to the hyperfine splitting 
of the high-lying states. Vacuum polarization effects can 
be neglected for states with £ > 2 to the order relevant 
for the current investigation. 

We split the calculation into a low-energy part, which 
contains Bethe logarithm type corrections (Sec. IHI[) . and 
a high-energy part, which can be treated on the basis of 
electron form factors fSec. IIVI) . For the low-energy part. 



we find that the net result can be expressed as the sum of 
corrections due to the hyperfine Hamiltonian, due to the 
energy correction, due to the wave function correction, 
and due to the hyperfine modification of the electron's 
transition current. For the high-energy part, we find a 
sum of two terms, one of which is due to a second-order 
effect involving the Pauli form factor correction to the 
Coulomb field, and the second of which is an anomalous 
magnetic moment correction to the hyperfine splitting, 
evaluated on relativistic wave functions. The first cor- 
rection can be split into two terms, which involve/do 
not involve mixing of the upper and lower components 
of the Dirac wave function, respectively. Quite surpris- 
ingly, the high-energy contribution can be expressed in 
closed analytic form, valid for an arbitrary excited state 
[sec Eq. ([55)) ]. For the Bethe logarithm type corrections 
relevant to the low-energy part, a numerical approach is 
indispensable. 

Finally, as indicated in Tables HI— IHIl we also find re- 
sults for D, F, and G states which arc of general interest 
to high-precision spectroscopy. 
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